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The amplification of magnetic fields in a highly conducting fluid is studied numerically. During 
growth, the magnetic field is spatially intermittent: it does not uniformly fill the volume, but is 
concentrated in long thin folded structures. Contrary to a commonly held view, intermittency of 
the folded field does not increase indefinitely throughout the growth stage if diffusion is present. 
Instead, as we show, the probability-density function (PDF) of the field strength becomes self- 
similar. The normalized moments increase with magnetic Prandtl number in a powerlike fashion. 
We argue that the self-similarity is to be expected with a finite flow scale and system size. In the 
nonlinear saturated state, intermittency is reduced and the PDF is exponential. Parallels are noted 
with self-similar behavior recently observed for passive-scalar mixing and for map dynamos. 

PACS numbers: 47.27.Gs, 91.25.Cw, 47.27.Eq, 47.65. +a, 95.30.Qd 
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We consider the problem of magnetic-energy amplifica- 
tion by a homogeneous isotropic turbulence in a conduct- 
ing fluid. This effect is known as the small-scale turbulent 
dynamo, and it is important to the dynamics of magnetic 
fields in astrophysical objects. It is, in fact, a generic 
property of random (in time and/or space) flows that 
they can amplify magnetic fluctuations at scales smaller 
than the scale of the flow itself. The amplification is a net 
result of stretching of the field lines by the local random 
strain associated with the flow (lj,|j,|3|- The fields gener- 
ated by this mechanism have a characteristic structure: 
they concentrate in flux folds containing antiparallel field 
lines that reverse direction at the resistive scale and re- 
main straight up to the flow scale 0, . ^ is because 
of the direction reversals that the magnetic energy in the 
wave-number space is predominantly at the resistive scale 
0] • The growing folded fields also concentrate in small 
parts of the system's volume: a phenomenon called inter- 
mittency. These properties of the small-scale dynamo are 
most pronounced for systems where the fluid viscosity is 
much larger than magnetic diffusivity (magnetic Prandtl 
number Pr m = v/rj 1), i.e., where magnetic fields re- 
verse direction at scales much smaller than the viscous 
cutoff. This regime is realized in many astrophysical plas- 
mas: examples are warm interstellar medium, intraclus- 
ter and intergalactic plasmas. In this Letter, we study the 
volume-filling properties of the dynamo-generated fields, 
i.e., the distribution of the field strength. 

Consider the equations of incompressible MHD: 



— u = z/Au- Vp + B- VB + f, V-u = 0, (1) 

at 



— B = B Vu + riAB, 

at 



(2) 



where p — const is density. Turbulence is excited by 
the forcing f . The spatial scale of the forcing is usually 
much larger than the diffusion scales of the fields. In 
our simulations, we solve Eqs. by a pseudospectral 

method. The forcing f is chosen to be random, white (<5- 
correlated) in time, and restricted to k/2ir = 1,2. The 
code units are based on box size 1 and mean injected 
power (u • f ) = 1. 

If diffusion is ignored (rj = 0), Eq. J5J) has the following 
formal solution in the comoving frame: 



]n[B(t)/B ] = / dt'bh : Vu(t'), 



where d/dt = dt +u- V. The pressure p and the magnetic 
field B are normalized by p and (^np) 1 ^ 2 , respectively, 



(3) 



where b = B/B. The Central Limit Theorem suggests 
that B should have a lognormal PDF with both the mean 
(In B) oc t and dispersion D cx t [cf. Eq. . This implies 
that (B n ) cx exp(7„i), where the growth rates depend 
quadratically on the order: 7„ oc n(n+3) 0,0.0]. Thus, 
in the diffusion-free case, the intermittency of the field- 
strength distribution increases in time in the sense that 
the kurtosis (B 4 )/(B 2 ) 2 and all other normalized mo- 
ments (B mn ) I (B m ) (B n ) grow exponentially. In space, 
intermittency means that the growing fields do not uni- 
formly fill the volume (compared to, e.g., a Gaussian field 
with the same energy). Assuming the equivalence of en- 
semble and volume averages, we may roughly interpret 
(B A ) / (B 2 ) 2 as an inverse volume-filling fraction. 

In the limit of small but finite rj, the small-scale dy- 
namo can still operate, but analytical description is much 
harder than for the diffusion-free case. The traditional 
approach that makes Eq. (J2J solvable in the diffusive 
regime is to consider a model random velocity field that 
is linear in space and white in time 0. The mo- 
ments of B can then be related to the distribution of 
the finite-time Lyapunov exponents associated with the 
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FIG. 1: Evolution of normalized moments of the magnetic- 
field strength. The (square root of) kurtosis of the velocity 
field is given for comparison. The kinematic (growth) and 
nonlinear (saturation) stages are demarcated in the plot. 



velocity-gradient matrix Vu, which is a function of time 
only The result is that the growth rates of (B n ) 

still increase with n as 7™ ~ n 2 , i.e, intermittency contin- 
ues to grow with time Q • It has so far been an accepted 
view that this model adequately describes the turbulent 
dynamo in the large-Pr m limit. In fact, the picture of 
increasing intermittency is not borne out by numerical 
experiments. 

In order to understand why, it is crucial to appreciate 
that the results obtained in the linear-velocity model ap- 
ply only as long as magnetic fluctuations are "unaware" 
of the fmitcness of the flow scale (system size). The most 
intuitive, albeit nonrigorous, argument as to why finite- 
scale effects should be important is as follows. 

The fields everywhere are stretched exponentially, but 
with fluctuating stretching rates, so any occasional dif- 
ference in field strength between different substructures 
tends to be amplified exponentially. Intermittency can 
grow with time if the system is infinite because for each 
moment (B n ), an ever smaller set of substructures can 
always be found in which the field has exponentially out- 
grown the rest of the system and which, therefore, dom- 
inantly contribute to (B n ). By contrast, in a finite sys- 
tem, only a finite number of exponentially growing sub- 
structures can exist, so the contribution to all moments 
must eventually come from the same fastest-growing one. 
The statistics of B should, therefore, be self-similar, with 
(B n ) growing at rates proportional to n, not n 2 , and all 
normalized moments (B mn )/(B m ) (B n ) saturating [2fj| . 



This is exactly what happens in our simulations. Af- 
ter initial diffusion-free growth, the normalized moments 
saturate (Fig.0. The PDF of the magnetic-field strength 
becomes self-similar: namely, the PDFs of B/B lms (here 
B lms = {B 2 ) 1 / 2 ) collapse onto a single stationary profile 
throughout the kinematic stage of the dynamo (Fig.|2J). 
The large- B tail of the PDF of B/B ims is reasonably well 
fitted by a lognormal distribution. Namely, suppose the 
PDF ofz= In B is 



P z {z) = (7nD)~ 1/2 exp[-(2- QnB)) 2 /D]. 



(4) 



D 



Then (B n ) oc exp [(\nB)n + Dn 2 /4] , so 
ln(( J B 4 ) 1 / 2 /(B 2 )). In the diffusive regime, D = const 
and the field-strength statistics become self-similar: the 
PDF of C = hi(-B/-B r ms) is stationary: 



P C (C) = frDrWexpl-k + Dpy/D]. 



(5) 



The lognormal fit in Fig. |3] is obtained by calculating D 
from the numerical data and comparing the profile (0 
with the numerically calculated PDF. The fit is qualita- 
tive but decent, considering the simplicity of the cho- 
sen profile JSJ, large statistical errors in determining 
(B 4 )/(B 2 ) 2 (Fig. |3J|, and dealiasing-induced numerical 
errors in resolving the field structure |6( . The PDF at low 
values of B appears to be power like (Fig. Eb), but there 
may be an unresolved lognormal tail at even smaller B. 
Note that the large- B tail describes the straight segments 
of the folds, while the small- B tail gives the field-strength 
distribution for the weak fields in the bends 0. 

The plots in Figs. are for a 256 3 simulation of 
Eqs. (PEl with v = 5 x 1CT 2 and Pr m = v/r) = 500. The 
Reynolds number for this run is Re = (it 2 ) 1 / 2 /vko ~ 2, 
(here fco = 2ir is the box wavenumber). Thus, the veloc- 
ity field, while random, is smooth in space. This is the 
so-called viscosity-dominated regime, which is the only 
physical setting in which large Pr m can be resolved nu- 
merically. Interestingly, results for larger Re ("real tur- 
bulence") and Pr m > 1 are very similar, especially in 
the kinematic regime. The reason for this is that the 
small-scale dynamo is always driven by the fastest eddies 
- the viscous-scale ones |s| , — and essentially the same 
field-stretching mechanism applies in both synthetic one- 
scale dynamos 0, S 0, 0] and in turbulent systems 
with Pr m > 1. In runs with Pr m = 10, Re ~ 100, and 
even with Pr m = 1, Re ~ 450 (not shown), we have 
found behavior analogous to that described above: satu- 
ration of normalized moments in the kinematic diffusive 
regime and self-similar field-strength PDFs fairly well fit- 
ted by the lognormal profile (J5J. In fact, the lognormal 
fit worked even better in these cases, which are some- 
what less violently fluctuating, because velocity is not as 
strongly coupled to the forcing. More detailed compar- 
isons will be reported in a future paper Q. 

The hypothetical lognormal PDF becomes self- 
similar only if its dispersion D does not depend on time. 
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FIG. 2: (a) Evolution of the PDF of the magnetic-field strength, (b) The PDFs of B / B lms in the kinematic regime: collapse 
onto a self-similar profile. The PDF in the saturated state is also shown. 



In contrast, in the diffusion-free regime, D ~ jt, where 
7 ~ Vu is the stretching rate (the turnover rate of the 
viscous-scale eddies). In the case of 77 > and linear ve- 
locity field, the formula for (B n ) derived by Chertkov et 
al. [9| is also consistent with a lognormal distribution for 
which D ~ jt. Both results are only valid transiently, 
during the time that it takes magnetic fluctuations to 
reach the resistive scale (in the former case) or the system 

(flow) scale (in the latter case). Since the scale separation 
1/2 

is ~ Pr^ and and the spreading over scales proceeds ex- 
ponentially fast at the rate ~ 7 (e.g. ||), the time during 
which intermittency increases is t* ~ 7 _1 InPr^ 2 . Physi- 
cally, this is the time necessary to form a typical fold with 
length of the order of the flow scale and field reversals at 
the resistive scale — starting either from a flow-scale or 
a resistive-scale fluctuation. We might conjecture that 

this time determines the magnitude of the dispersion in 

1/2 

the self-similar regime: D ~ t** ~ InPC , which im- 
plies that the kurtosis (B 4 )/(B 2 ) 2 = exp(2D) increases 
with Pr m in a powerlike fashion. The specific power law 
depends on prefactors that may be non-universal. The 
same holds for other normalized moments. Figure^is an 
attempt to test this hypothesis for a sequence of simula- 
tions with increasing Pr m . The fluctuations in the kine- 
matic regime are very large (see Fig.^l, so the error bars 
are too wide to allow us to claim definite confirmation 
of the powerlike behavior, but our results are consistent 
with a (£ 4 )/(B 2 ) 2 - Pr°; 3 scaling HJ. 

We emphasize that the self-similarity reported here is 
statistical, not exact. Namely, it does not imply that 
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FIG. 3: The PDF of B/B ruls averaged over the kinematic 
diffusive stage and the lognormal profile with the same D. 
Also given is the PDF in the saturated state. 



the magnetic field is simply a growing eigenmode of the 
induction equation J5J. Such an eigenmode does exist 
for some finite-scale non-random (and time-independent) 
flows and maps, owing essentially to the fact that the dif- 
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FIG. 4: The kurtosis {B 4 )/(B 2 ) 2 during the diffusive kine- 
matic and the nonlinear stages of a sequence of runs with 
v = 5 x 1CT 2 (low Re) and 25 < Pr m < 2500. All averages 
are over 20 time units (~ 20 turnover times). 



fusion operator has a discrete spectrum in a finite domain 
Q . The self-similar PDF we have found here is a natural 
counterpart of this eigenmode dynamo for random flows. 

Note that for the problem of passive-scalar decay, 
self-similar behavior was also found in certain two- 
dimensional maps (the so-called "strange mode" 0, 0, 
ITS llijp. two-dimensional non-smooth inverse-cascading 
turbulence ^J, and even in scalar-mixing experi- 
ments [l^- In map dynamos studied by Ott et al. [SIf^. 
moments of B also grew at rates cx n (these authors re- 
lated such behavior to the flux-cancellation property of 
the field, i.e, to the folded structure). We expect that 
self-similar evolution is a fundamental property of pas- 
sive advection of scalar and vector fields by finite-scale 
flows. 

A detailed discussion of the nonlinear saturation of the 
dynamo and of the magnetic-field intermittency in the 
saturated state is beyond the scope of this Letter. We 
limit ourselves to mentioning that the nonlincarity leads 
to a reduction of intermittency due to tighter packing of 
the system domain by the saturated fields (Fig. [TJ. It 
is not a surprising result: nonlinear back reaction im- 
poses an upper bound on the field growth, and, once the 
strongest fields in the dominant substructure saturate, 
the weaker ones elsewhere have an opportunity to catch 
up. The PDF of the saturated field turns out to be expo- 
nential (Fig. El cf. Olli). Accordingly, in the nonlinear 
regime, the kurtosis does not depend on Pr m (Fig. |3J. 
Furthermore, while B Tms is not Pr m -independent for fi- 



nite values of Pr m (though B rms does tend to a constant 
value ~ u rms as Pr m — >• oo), the PDF of B/B Tms turns 
out to be the same for all Pr m . Further results on the 
nonlinear dynamo will be reported elsewhere Q . 
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Note that intermittency is often associated with the pres- 
ence of sparse field structures that have disparate spatial 
dimensions ("coherent structures"). In fluid turbulence, 
these are the vortex filaments. For the dynamo-generated 
magnetic fields, the coherent structures are the folds, 
for which the disparate dimensions are their length and 
the field-reversal scale (~ resistive scale). In the linear- 
velocity model, the folds are elongated indefinitely. In 
reality, the length of the folds cannot be larger than the 
scale of the flow because of the bending jjj. The scale 
separation in the folds is, therefore, bounded from above 
by ~Pri( 2 . 

A note of caution is in order: at current resolutions, our 
lognormal fit is not the unique possibility consistent with 
numerical results: stretched-exponential and even steep- 
power-tail fits of comparable quality can be achieved. 
Note, however, that a stretched exponential would not 
be consistent with Pr m -dependent normalized moments. 



